International Journal of Engineering and Applied Sciences (IJEAS) 

ISSN: 2394-3661, Volume-1, Issue-3, December 2014 


A Theoretical Model for the Design of a Labyrinth 

Seal Gland 

E. Saber, Khaled M. Abdou 


Abstract — The analysis of a rectangular labyrinth seals 
having several configurations is carried out in this study. 
The laminar, compressible and axisymmetric flow in a 
rectangular labyrinth seal is treated and the seal 
performance is obtained. The analysis is based on the 
assumptions that adiabatic conditions exist and that the 
ratio of the seal height to the shaft radius is small as 
always the case in turbomachines. The theoretical work 
is capable of handling single cavity in order to compute 
the pressure ratio, loss in stagnation pressure, change in 
density, temperature ratio and the variation of inlet to 
exit Mach number across the cavity. These quantities are 
obtained at several values of leakage flow Mach number 
to be used in estimating the required number of cavities 
to get a minimum allowable leakage rate. For a subsonic 
leakage flow that is not in choked condition, it is 
concluded in this work that the performance relations of 
a single cavity can be used to design a multi -cavity seal to 
meet an allowable leakage rate. 

Index Terms — Labyrinth seals - Non-contact seal - 
Turbomachines - Compressors - Computational fluid 
mechanics - Mathematical models of fluid flow. 


I. INTRODUCTION 

Effective sealing in turbomachines is required in order to 
minimize the leakage of high-energy working fluid, thus 
reducing the adverse effects of seal leakage on overall 
engine performance. Despite the development of new 
techniques in sealing technology, labyrinth seals will 
continue to be the most widely used sealing element in high- 
speed turbomachines. Simplicity and compactness, 
reliability and temperature resistance are some of its major 
advantages. In general, increasing the fluid resistance 
through the labyrinth seal can reduce leakage. This can be 
accomplished either by reducing the clearance size, the 
development of a new geometrical seal configuration and the 
use of blowing through the seal gland. Many investigators 
tackled the leakage problem through labyrinth seals for 
turbomachinery. Some of the early design methods available 
in the literature are purely analytical [1-3] and others are 
analytical-experimental [4-7]. Benvenuti et. al. made a 
comparative study of these methods and they found that such 
treatments provided results with considerable differences 
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and contradictions. The first reliable prediction method was 
due to Stoff [8] who reported a theoretical and experimental 
study of incompressible turbulent flow through a labyrinth 
seal of the straight-through grooved shaft type. He was able 
to show how much the flow pattern in the labyrinth seal is 
affected by the ratio of the axial leakage velocity to the shaft 
peripheral speed. It was found by Rhode et. al. [9] that 
labyrinth seal solutions using the difference scheme of 
Stoff [8] may suffer from the so-called false diffusion 
numerical error depending on the flow field conditions. They 
however, presented an approach to overcome the difficulties 
encountered when using the scheme of Stoff. They also 
adopted for a typical cavity configuration a simple stair step 
approximation to the curved solid walls. El-Gamal et. al. 
[10-12] were able to examine the pressure drop and flow 
pattern for a number of geometrical configurations of seals. 
Their results showed how important the effect of shaft 
rotation is on the performance of the seal. Yizhang and Feng 
[13] applied the continuity equation, energy equation, 
Bernoulli's equation and the equation of state to the vortex 
cavity to obtain a closed set of equations. They solved these 
equations simultaneously and numerically to determine the 
velocity, pressure and temperature in each vortex cavity. 
Computational fluid dynamics (CFD) is increasingly used in 
analyzing flows inside labyrinth seals. An advantage of 
using CFD is its capability to analyze a large number of 
design configuration parameters in a relatively short period 
of time. Therefore, with the development of commercial 
codes, the use of CFD analysis has been increasing rapidly 
in recent years as suggested in the references [14-18]. Tong 
and Kyu [19] analyzed the influence of configuration and 
clearance on the leakage behavior of straight and steeped 
labyrinth seals with various clearances. Both computational 
fluid dynamic (CFD) and an analytical tool were used to 
predict the leakage flow through the straight and steeped 
labyrinth seals. They were compared the predicted results 
with experimental data [14]. The CFD predicted sufficiently 
well the leakage flows of different seal configurations and 
for various clearances, especially for the stepped seal. The 
analytical model provided acceptable prediction in the case 
of straight seal, but did not provide satisfactory predictions 
for the stepped seal. They stated that a straight seal with 
more teeth seems to be a suitable design. Zhao et al [20] 
studied the effects that cavity dimensions, numbers, and 
locations have on the leakage flow in straight-through 
labyrinth seals. The Reynolds averaged Navier-Stokes 
equations along with a standard k — CD turbulence model 
were solved and the flow is assumed to be axisymmetric. 
They found that the cavity dimensions and cavity number 
have significant effects on the leakage as well as on the flow 
pattern in the seal. The aim of the present work is to 
introduce a theoretical model to be used in estimating the 
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required number of cavities to get a minimum allowable 
leakage rate. 

II. ANALYSIS 

The geometrical configuration of the labyrinth seal selected 
and the coordinate system used in the analysis are shown in 
Fig. 1. In the absence of entry swirl, the dimensionless 
equations governing the steady laminar compressible 
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u , v and w are the velocity components in dimensionless 
form in the 6,y and z, directions respectively. 

p,p ,Ju and A are the dimensionless pressure, density, 
viscosity and dissipation function respectively. 
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R e , E, M and P r are the Reynolds number, Eckert number, 

Mach number and Prandtl number respectively. Also, we can 
express the Eckert number in term of Mach number as, 

E = M 2 (y g -l). 

Let, v =o( m), w = o( 1 ), u = o( 1 ), p = 1, 
p =o( 1 ), y -o( 1 ) and z -o( 1 ) ( 8 ) 

Assuming that the pressure is constant across the seal height 
and applying the order of magnitude of (8) into equations (1) 
, (2) , (3) , and (4), then 


P2 ~ 
o u 

-2 

dy 


2 ou — _ — du 

+ y — r = ypR e w -= 

dz 


-2 

dz 


e -l=o 

dy 


dw d p 

y p w — = —y ——r 
dz dz 


(9) 


( 10 ) 


R. 


d z w 4 2 d z w 

'■—y 

3 


-2 

dy 


-2 

dz 


( 11 ) 


— dw ~dp y ^ v 

p + w = 0 (12) 

dz dz 

Assuming adiabatic flow condition and neglecting the 
temperature variation across the labyrinth height, T =T( z), 
p = p(z) and equations (5) and (7) can be written as 
follows, 
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Integrating equations (13) and (14) with respect to y , we 
get 
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Using equation (6) and differentiate with respect to Z , the 
variation of the pressure in axial direction can be obtained as, 
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Integrating the continuity equation (12) with respect to y and 
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using I wdy - G(z) , we get the continuity equation in the 


following form, 
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Here, A is a constant to be calculated from: A = G (0) 
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The grooved shaft is the only case treated here and presented 
as an example of a class of rectangular labyrinth seals. The 
presence of a clearance in the labyrinth seal causes an axial 
leakage of the fluid, this means that the axial leakage velocity 
component must be taken into consideration. It is assumed 
that the expansion through the clearance is isentropic and that 
no slip condition at the walls. The boundary conditions for the 
grooved shaft labyrinth shape selected here are, 
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III. Leakage Model 

The present model may be employed to determine the total 
number of labyrinth cavities for a wide variety of designs. It 
does, however, require results from single-cavity 
computations. Hence, the prediction of leakage for a new 
design involves two step procedure. First, the single-cavity 
solution is used to determine the inlet to outlet cavity 
pressure ratio, temperature ratio and Mach number for the 
given seal configuration. This is done at several cavity inlet 
Mach numbers and at a fixed values of inlet stagnation 
pressure and stagnation temperature. Second, the inlet-to- 
outlet cavity pressure ratio, temperature ratio and 
subsequently the cavity outlet static pressure and static 
temperature are evaluated from the polynomial expressions 
at the known inlet Mach number. The cavity outlet density is 
evaluated from the ideal gas law. The cavity outlet velocity 
and subsequently the corresponding Mach number is 


calculated from the mass conservation equation. Since the 
outlet conditions for a given cavity are the inlet conditions 
for the next downstream, the same sequence is repeated for 
each cavity. At last cavity, the exit pressure must be less than 
or equal downstream pressure to estimate the total number of 
cavities performed the seal gland. The leakage model was 
developed as a simple computer program to which the user 
specifies the desired mass leakage rate, the upstream 
stagnation conditions, the downstream pressure, shaft 
diameter, shaft speed and gland diameter. 


IV. RESULTS AND DISCUSSION 


The governing equations (9), (15), (16) and (19) with the 
boundary conditions for the flow in grooved shaft 
rectangular labyrinth seals (20) are suited to numerical 
solution by the finite difference technique using the Gauss- 
Seidel under-relaxation iterative method. The domain of 
integration is divided into 40 x 40 square grids each of 
mesh size h = 0.025. The results are obtained for a single 
grooved shaft cavity and for different values of aspect ratio 
y . It is found that the variations in results for y = 0( 1 ) are 

less than 1% and may be neglected. Thus the results for the 
steady flow in grooved shaft labyrinth seals having y -1 

only are presented. 

The results thus obtained are used to study the variation of 
various parameters and conditions of flow at the inlet and 
exit of a single straight through grooved shaft labyrinth 
cavity. The performance curves and relations of a single 
cavity are obtained for different values of the speed of 
rotation of the shaft. Solving equation (16) gives the 
temperature distribution inside the cavity, thus the Mach 
number at exit can be calculated as 





Fig. 2 shows the variation of the exit Mach number M e with 
the inlet Mach number M t for different values of the 


rotational Mach number M r , M r 



For a 


certain value of the inlet Mach number, it is found that 
increasing the rotational Mach number leads to decreasing 
the value of the exit Mach number. 


P 

Fig. 3 presents the relation between the pressure ratio — 

Pi 

and the inlet Mach number for different values of 

rotational Mach number M r . The figure shows a large 

pressure drop at rotational Mach number M r =0.125, 

increasing the speed of the shaft the rotational Mach number 
increases and the pressure drop across the seal decreases. 
This means that increasing the speed of the shaft decreases 
the acceleration of the flow in the axial direction. 


In the design of a cavity of a labyrinth seal, the value of the 
downstream pressure is significant. The variation of the ratio 
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of the stagnation pressure at inlet to static pressure at exit 

p • 

— — with the inlet Mach number M t is shown in Fig. 4. 

Pe 

Also, Fig. 5 shows the variation of the static temperature 
T 

ratio — with the inlet Mach number M ; for different 

rri t 

* i 

values of rotational Mach number M r . It is found that 

increasing the speed of the shaft decreases the change in the 
static temperature . Also, the maximum inlet Mach number 
the maximum decrease in static temperature across the 
cavity. The leakage rate through a single cavity seal depends 
on the inlet conditions, which are stagnation pressure and 
stagnation temperature, p oi and T oi and the value of the 

exit static pressure p e . Introducing the mass (leakage) flow 

mJRT oi 

parameter as co D . T . = — - , we get a useful relation 

F O l O l m 4 

Poi ■ C 

between the leakage per unit area and stagnation temperature 
and stagnation pressure at inlet. Fig. 6 shows the variation of 

p 

leakage flow parameter c ° Poi T oi w hh pressure ratio — 

Pi 

across the seal cavity for various rotational Mach numbers. 
The results show that increasing the rotational Mach number 
in general leads to an increase in the leakage rate. 


For multi cavity seal, p ou and T ou are the upstream 

stagnation pressure and upstream stagnation temperature 
respectively. Also, p u and p d are the upstream and 


downstream static pressures at inlet and exit of a multi cavity 
seal respectively. Defining the leakage rate parameter 

rhjRT ou 

through a multi cavity seal as co v T - — - , Fig. 7 

F O 14 O 14 4 

Pou Ac 

shows the relation between the leakage rate parameter 


p 

co n r and the pressure ratio — — for different number of 

you^ou A „ 

Pd 

cavities that formed a multi cavity seal. The results show that 
the leakage rate decreases as the number of cavity increases. 
Fig. 8 shows the effect of number of cavities on the leakage 
rate parameter for different values of pressure ratio of 
upstream stagnation pressure to downstream pressure. It is 
seen that the larger the number of cavities the smaller the 

p 

leakage rate parameter at any pressure ratio — 1 — . However 

Pd 

there is an asymptotic value of number of cavities at which 
leakage rate parameter c °p ol j ou becomes very small. 


Gland diameter =88 mm 
Allowable leakage rate =0.001 kg/s 

It is required to make a preliminary design to estimate the 
required number of cavities of the labyrinth seal gland. 


Design procedure: 

1- Inlet Conditions of the first cavity: 

Assuming isentropic expansion through the inlet clearance 
and using the mass conservation equation we obtain the inlet 
Mach number, inlet static pressure and inlet static 
temperature. 

2- Polynomial expressions 

The resulting inlet to outlet single cavity pressure ratio, 
temperature ratio and exit Mach number are inserted into the 
leakage model using polynomial representation versus inlet 
Mach number. The polynomial expression for the single 
cavity performance curves may be written as, 
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Where, 

The coefficients of inlet/exit Mach number relation are, 
a Q =0.000636163, aj = 1.12124, a 2 =0.122316 

The coefficients of inlet/exit pressure ratio relation are, 
b Q = 1.00027, bj = -2.60013, b 2 = 19.7262, b 3 = -78.5611 

b 4 = 175.898, b 5 = -224.228, b 6 = 151.756, b 7 = -42.3422 

The coefficients of inlet/exit temperature ratio relation are, 
^ = 0.999979, ej = -0.00275745, e 2 = -0.0491327, 

and e 3 - -0.0375619 


3- Exit conditions of the first cavity 

Using polynomial expressions, the exit Mach number of the 


first cavity Mj , the pressure ratio 


V 

Pi J 


and temperature 


ratio 


{ rj, 

L l 


T 

V A 7 


can be calculated. Also, the exit static pressure 


( pi ) and static temperature ( 7} ) may be estimated. 


I. CASE STUDY 

Consider a labyrinth seal of air copmpressor with the 
following conditions: 

Upstream stagnation pressure = 140 kPa 

P 

Pressure ratio — — =1.4 

Pd 

Up stream stagnation temperature = 400 K 

a) Shaft diameter = 80 mm 
Shaft speed = 200 rps 


4- Estimation of the total number of cavities 

(1) Total number of cavities =11 

Estimated mass flow rate= 0.00081 12 kg/s 
The pressure distribution along the multi cavity seal 
obtained and presented in Fig. 9. Also, Fig. 9 shows the 
pressure distribution along the labyrinth seal for different 
values of allawable leakage rate, it is found that increasing 
the number of cavities decreases the leakage rate 
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V. CONCLUSIONS 


A simple approach to the laminar compressible flow 
problem through labyrinth seals of rectangular shapes has 
been developed. The theoretical model enabled the effect of 
shaft rotation on the performance of the seal to be included. 
It is found that increasing the speed of the shaft decreases 
the pressure drop across the seal. For a single cavity and at 
several leakage flow Mach number, the pressure ratio, 
temperature ratio and the variation of inlet to exit Mach 
number across the seal can be obtained. These quantities are 
used to estimate the required number of cavities to get a 
minimum allowable leakage rate. For a subsonic leakage 
flow that is not in choked condition, it is concluded that the 
performance relations of a single cavity can be used to 
design multi-cavity seal to meet an allowable leakage rate. 
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Fig. 1 . Labyrinth seal coordinate system. 



Fig. 2. Variation of seal inlet to exit Mach number. 



Fig. 3. Variation of static pressure ratio across 
the labyrinth cavity with inlet Mach number. 
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Fig. 4 Pressure ratio p oi / p e against 
the inlet Mach number 
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Fig. 5. Variation of temperature ratio — 
with inlet Mach number. 



Fig. 6. Leakage rate parameter &> p t c 


against pressure ratio — 

Pi 
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Fig. 9. Pressure distribution along the multi-cavity labyrinth 

seal. 
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